"Light Sail" Acceleration Revisited 
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The dynamics of the acceleration of ultrathin foil targets by the radiation pressure of superin- 
tense, circularly polarized laser pulses is investigated by analytical modeling and particle-in-cell 
simulations. By addressing self-induced transparency and charge separation effects, it is shown that 
for "optimal" values of the foil thickness only a thin layer at the rear side is accelerated by radia- 
tion pressure. The simple "Light Sail" model gives a good estimate of the energy per nucleon, but 
overstimates the conversion efficiency of laser energy into monoenergetic ions. 
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Radiation Pressure Acceleration (RPA) of ultrathin 
solid targets by superintense laser pulses has been pro- 
posed as a promising way to accelerate large numbers of 
ions up to "relativistic" energies, i.e. in the GeV/nucleon 
range 0, 0, i, II H, @, 0, B H ■ The simplest model of this 
acceleration regime is that of a "perfect" (i.e. totally re- 
flecting) plane mirror boosted by a light wave at perpen- 
dicular incidence [Io| . which is also known as the "Light 
Sail" (LS) model. The LS model predicts the efficiency 
77, defined as the ratio between the mechanical energy of 
the mirror over the electromagnetic energy of the light 
wave pulse, to be given by 



77 = 2/3/(1 + 0), (3 = V/c, 



(1) 



where V is the mirror velocity; hence, RPA becomes 
more and more efficient (77 — > 1) as /3 — > 1. Heuristi- 
cally, Eq. |T]) can be explained by the conservation of the 
number of "photons" N of the light wave reflected by 
the moving mirror in a small time interval: each pho- 
ton has energy fiuj, thus the total energy of the incident 
and reflected pulses are given by Nhu> and NfiLo r , where 
u> r = ui(l — /3)/(l +/3) due to the Doppler effect, and the 
energy transfered to the mirror is given by their difference 
[2/3/(1 + 0)]NHu. 

The predictions of the LS model are very appealing 
for applications, but one may wonder to what extent this 
picture is appropriate to describe the acceleration of a 
solid target by a superintense laser pulse. In the present 
paper, we revisit the LS model with the help of simple 
modeling and particle-in-cell (PIC) simulations. We ad- 
dress issues outside the model itself, such as the effects of 
nonlinear reflectivity and charge depletion, and on this 
basis we explain a few features observed in simulations. 
Our main result is that the LS model is accurate in pre- 
dicting the ion energy but overstimates the corresponding 
conversion efficiency, i.e. the fraction of the laser pulse 
energy transferred into quasi-monoenergetic ions, due to 
the fact that only a layer of the foil at its rear side is 
accelerated by RPA. 

Our analysis is confined to a one-dimensional (ID) 
approach for the sake of simplicity and because multi- 
dimensional simulations showed that a "quasi- ID" ge- 
ometry has to be preserved in the acceleration stage 



(by using flat-top intensity profiles) to avoid early pulse 
transmission due to the expansion of the foil in the ra- 
dial direction Circularly polarized pulses are used 



to reduce electron heating 
by several groups 

E3. 



12| . an approach followed 



Dy several groups for efficient acceleration of thin foils 
[1, 01, 0, 0, @, EH • We do not consider intensities high 
enough that ions become relativistic within the first laser 
cycle; this condition may affect the early stage of charge 
depletion (e.g. by narrowing the temporal scale separa- 
tion between ions and electrons), and lead to different 
estimates [H, 0]. 

The LS model is based on the following equation of 
motion for the foil 



(2) 



where 7 = 1/^/1 — /3 2 , dX/dt — V, I is the light wave 
intensity, p and £ are the mass density and thickness of 
the foil, R(u>') is the reflectivity in the rest frame of the 
foil, and u>' = u> v / (l-/3)/(l+/3). For suitable expres- 
sions of R(oj'), the final velocity (3 t can be obtained from 
Eq.@ as a function of the pulse fluence F = J Idt. For 
R = 1, one obtains 
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(l + £) 2 + l' ~ pic 2 " Am 

In the last equality we wrote the fluence in dimension- 
less units as a^r, where ao = y^I/m e c 3 n c is the dimcn- 
sionless pulse amplitude and r is the pulse duration in 
units of the laser period, and introduced the parameter 
£ = Tt(no/n c )(£/\) which characterizes the optical prop- 
erties of a sub-wavelength plasma foil 13]. In these equa- 
tions, no is the initial electron density, n c = 7rm e c 2 /e 2 A 2 
is the cut-off density, and A is the laser wavelength. 
In practical units, n c — 1.1 x 10 21 cm~ 3 [A//zm] -2 and 
a = (0.85/V2)(/A 2 /10 18 W cm^pm 2 ) 1 / 2 for a circu- 
larly polarized laser pulse. Using Eq.Q it is found that 
with a 1 ps = 10~ 12 s, 1 PW = 10 12 W laser pulse and a 
10 nm target of 1 g cm' 3 density, ~ 1 GeV per nucleon 
may be obtained. As the LS model assumes the target 
to be a perfect mirror (i.e. rigid and totally reflecting), 
it implies that all the ions are accelerated to the same 
velocity and the spectrum is perfectly monoenergetic. 
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FIG. 1: (Color online) Parametric study of the ion energy 
spectra vs. laser amplitude ao and foil thickness I. The 
contours of log 10 fi(E) are shown, with fi(E) the energy 
per nucleon distribution normalized to unity. For all runs, 
no = 250n c , Z/A = 1/2, r = 9. The dashed line shows the 
prediction of the LS model for the ion energy. The dotted 
horizontal line marks £ opt given by the f = ao condition. 



FigQ] shows a parametric study of the ion spectrum vs. 
I and ao from PIC simulations. For all runs, uq = 250n c , 
Z/A = 1/2 and the pulse has a flat -top envelope with 1 
cycle rise and fall times and 8 cycles plateau. For each 
value of ao and for £ less than a threshold value £ opt we 
observe a narrow spectral peak, whose energy increases 
with decreasing £ and is in very good agreement with the 
predictions of the LS model, assuming R = 1. A typical 
lineout of the spectrum is shown in Fig[2]a). For £ > 
£ opt , the peak disappears and a thermal-like spectrum 
is observed. This is correlated with an almost complete 
expulsion of the electrons from the foil in the forward 
direction at the beginning of the interaction, leading to 
a Coulomb explosion of the ions. 

The results of FigQ] show that the LS model is useful 
for quantitative predictions of the ion energy, but also 
suggest several questions of interest both for the basic 
physics of RPA and its applications. How is £ opt deter- 
mined? Does the reflectivity of the foil and relativistic 
effects on the latter play a role? As the radiation pressure 
tends to separate electrons from ions, does the foil re- 
main neutral before and/or after the acceleration stage? 
Moreover, as shown shown in Fig[2]b), the "monoener- 
getic" peak contains just a fraction of the total number 
of ions, and such fraction depends on £ and a . This 
is different from the assumption of the LS model, which 
assumes all the ions in the foil to move coherently with 
the foil, and may sound surprising, since the peak energy 
is in agreement with the LS formula where the whole 
mass of the foil, including low-energy ions out of spectral 
peak, is used. In the following we provide answers to the 
questions above by discussing effects not included in the 
simplest LS model, i.e. beyond the description of the foil 
as a perfect, rigid mirror. 

First we discuss effects related to the reflectivity R of 
the plasma foil. For very high intensities, electrons oscil- 
late with relativistic momenta in the laser field, leading 
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FIG. 2: (Color online) a): Ion energy spectra (in energy per 
nucleon) from a simulation with ao = 30 and a I — 0.04A 
thick foil of a single ion species with Z/A — 1/2 (top) and 
one with the same parameters but where ions in a thin sur- 
face layer (0.01A) at the rear side are replaced by protons 
(bottom). Fig[T] b) Fraction of ions contained in the spectral 
peak vs. the target thickness I for three values of ao = 10 
(black, crosses), 30 (blue, triangles) and 50 (red, squares). 
The dashed lines correspond to Eq.© for F. All other pa- 
rameters for both a) and b) are the same as in Fig[T] c) Ap- 
proximate profiles of ion (m, green) and electron (n e , blue) 
densities and of the electrostatic field (E x , red) in the early 
stage of the interaction, before ions move. 



to a nonlinear dependence of R upon ao. An explicit ex- 
pression can be found analytically by using the model of 
a delta-like "thin foil" [13[, i.e. a plasma slab located 
at x — with electron density n e (x) = na£8(x). The 
expression obtained for R in the rest frame of the foil is 
very well approximated by 



R 



C 2 /(l + C 2 ) («o< y/T+e) 

C 2 /a 2 (a > VT+e) 



(4) 



A threshold for self-induced transparency of the foil may 
thus be defined as ao = yl + ( 2 — ( when ( 3> 1, i.e. 
in most cases of interest. According to Eq. ([!]), the total 
radiation pressure P lad on the target 



P rad = 2RI/c = 2m e c 3 n c alR 



(5) 



becomes independent upon ao for ao > C Thus, the 
maximum radiation pressure is obtained for ao ^ Ci an< i 
in this condition typically R ~ 1 for solid densities. This 
suggests that the optimal thickness £ opt is determined 
by the condition ao — Ci i n good agreement with the 
simulation results in FigQ] and as also found by other 
studies dEl. 

The nonlinear reflectivity of the thin foil is determined 
by the transverse motion of electrons (in the foil plane). 
However, for a thin but "real" target the radiation pres- 
sure tends to push electrons also in longitudinal direction, 
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and may remove them from the foil. Let us compare P rad 
with the electrostatic pressure P BS that would be gener- 
ated if all electrons would be removed from the foil. The 
condition 

P., d > P s = 2ir(en £) 2 (6) 

corresponds to the threshold for the removal of all elec- 
trons from the foil. However, when Eq.(j4]) is used for 
«o < \A + C 2 > Eq-© reduces to a > \J\ + C 2 , while for 
a o > yl + C 2 we fi- n< i that P rad = P OB holds. It is thus 
possible to produce a density distribution where all elec- 
trons pile at the rear surface of the foil. In fact, if ao ^ C 
and R ~ 1, the laser pulse compresses the electron layer 
while keeping R constant since the product n e £ does not 
change during the compression; at the same time almost 
no electrons are ejected from the rear side because the 
ponderomotive force vanishes there if R ~ 1, and so does 
the electrostatic field: the qualitative profiles of the elec- 
tron density and of the electric field are shown in Figj2]c). 
Since for ao close to C the equilibrium between the elec- 
trostatic and radiation pressures occurs only when the 
depth d of the region of charge depletion is close to £, 
electrons are compressed in a very thin layer. The de- 
pletion depth d may be estimated from the equilibrium 
condition 

P OB = 2n(en d) 2 ~ 2I/c (7) 

when R ~ 1, which yields d ~ £(ao/C) ~ 1- It is worth 
to point out that these considerations are appropriate for 
a circularly polarized laser pulse; for linear polarization, 
all electrons may be expelled for a transient stage under 
the action of the J x B force whose peak value per unit 
surface exceeds 2RI/c due to its oscillating component. 
Complete expulsion of electrons for ao > C has been dis- 
cussed in Ref. [14j]. 

The snaphshots from a PIC simulation shown in Fig[3] 
for a case with ao = 30 and ( = 31.4 confirm the scenario 
outlined above. The electron density n e reaches values 
(out of scale in Fig[3]) up to tens of the initial density. A 
very high resolution Ax — A/2000 is used to resolve the 
density spike properly. For a laser pulse with flat-top 
envelope the density spiking at the rear side of the foil 
is particularly evident, but we verified that it occurs also 
for a "sin 2 " envelope. Similar features were observed also 
in Refs.[B|, Q, but not discussed in detail. 

The electron compression in a thin layer during the ini- 
tial "hole boring" stage has important consequences for 
the later acceleration stage. Let us refer to the approxi- 
mate field profiles in the initial stage, sketched in Fig[2]c), 
which were the basis of the model presented in Ref.[l2j. 
This model suggests that only the ions located initially 
in the electron compression layer (d < x < £) will be 
bunched and undergo RPA (via a "cyclic" acceleration 
as discussed in Refs.jH, 0, 0|) because for these ions only 
the electrostatic pressure balances the radiation pressure, 
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FIG. 3: (Color online) Snapshots from a ID PIC simulation 
of the interaction of a laser pulse with a thin plasma slab. 
The ion density m (green), the electron density n e (blue), the 
longitudinal electric field E x (red, dashed) and the pulse field 
amplitude El = \l Ey + El (red, dotted) are shown. The 
target left boundary is at x = where the pulse impinges 
at t — 0. Times are normalized to the laser period T, fields 
to m e uic/e, and densities of n c . The laser pulse has ampli- 
tude ao = 30 and the foil thickness is I — 0.04A. All other 
parameters are the same of Fig[T] 

while the ions in the electron depletion layer (0 < x < d) 
will be accelerated via Coulomb explosion, i.e. by their 
own space-charge field. This is exactly what is observed 
in the PIC simulations, both in the density profiles (see 
Figj3]at t = 2.2T) and in the ion spectra. This effect also 
explains how RPA with circularly polarized pulses may 
work also in double layer targets [15| , if the thickness of 
a thin layer on the rear side matches £ ef( = £ — d. Fig[5] 
shows ion spectra for the same simulation of Fig [3] and for 
a simulation with the same parameters, but where ions 
in a surface layer of 0.01A thickness have been replaced 
by protons. A fraction of heavier ions is also accelerated 
to the same energy per nucleon as the protons, a typical 
feature of RPA of a thin plasma foil. 

As an additional consequence of the piling up of elec- 
trons at the rear surface, the portion of the foil which is 
boosted by the laser pulse is negatively charged due to 
the excess of electrons. However, the simulations show 
that when the laser pulse is over the excess electrons de- 
tach from the foil and move in the backward direction, 
so that the accelerated layer is eventually neutral. This 
is important to avoid a later Coulomb explosion of the 
layer and to preserve a monoenergetic spectrum. During 
the acceleration, the longitudinal field at the surface of 
the accelerated layer is almost constant implying that the 
charge there contained is also constant. It is thus pos- 
sible to estimate the fraction F of accelerated ions from 
the initial equilibrium condition, Eq.([7j), as 

F~£ o!! /£~(l-a /C)- (8) 

The agreement with data in Fig|2]b) is qualitative, with 
large deviations as F becomes significantly smaller than 
one. As explained below, a lower bound on F is deter- 
mined by energy conservation. 

A simple argument of force balance also explains why 
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the energy of the spectral peak in Fig[2] is in very good 
agreement with the predictions of the LS model where 
the initial value I of the foil thickness is used, while only 
a layer of thickness £ ef( < I is accelerated via RPA. Let 
us refer again to the profiles of FigJ^l The equilibrium 
condition for electrons implies 
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= / eE 

Id 



e-d 



i p dx = -enoEol. 



(9) 



The electric field pushes ions in the compression layer 
d < x < £, exerting a total pressure 



Pr = I eE 

Id 



£-x\ , 21 [ e-d 

— d r idx= -{— 



(10) 



where we used Eq. §§§ and assumed R = 1. The equation 
of motion for the ion layer, in the early stage, can be thus 
written as 



which is trivially equivalent to 

= 

at c 



(11) 



(12) 



i.e. to the equation of motion one would write for the 
whole foil. The argument may be applied also when the 
layer is in motion leading to the same conclusion. Having 
the same (3{t) as the whole foil implies that the energy 
per nucleon and the efficiency (p} will be also the same, 
but the total kinetic energy will be lower for the thin 
layer. The rest of the absorbed energy is stored in the 
electrostatic field and as kinetic energy of the ions in 
the x < X(t) region. Let us consider for example the 
energy stored in the electrostatic field. At the time t, 
the field E x between the initial [x = 0) and the actual 
(x = X(t)) positions of the front surface of the foil is 
given approximately by E x = Eqx/ X(t), where Eq = 
Airenod, corresponding to an electrostatic energy per unit 
surface 



U m = U m (t) = 



which varies in time as 



X[t) El{x,t) 



8tt 



dx. 



(13) 



(14) 



Dividing (fT4f by the laser intensity we obtain the "con- 
version efficiency" into electrostatic energy ?y es 



1C 

I dt 



(15) 



If £ ~ ao and thus d ~ £, we would obtain ?7 C3 ~ 2/3 > r\ 
that is unphysical. Thus, the energy stored in the electro- 
static field also prevents the accelerated layer thickness 
to shrink to zero. 



In conclusion, we have revisited the "Light Sail" model 
of Radiation Pressure Acceleration of a thin plasma foil. 
The nonlinear reflectivity of the foil determines the "op- 
timal" condition ( ~ a , for which the energy in the RPA 
spectral peak is highest and in good agreement with the 
LS model formula where the total thickness (or the total 
mass) of the foil enters as a parameter. However, not 
all the foil is accelerated, but only a thin layer at the 
rear side of thickness £ eff < i\ the apparent paradox is 
solved by observing that, to keep electrons in a mechan- 
ical quasi-equilibrium, the electrostatic pressure pushing 
ions in the accelerated layer is £ ett /£ times the radiation 
pressure on electrons, so that the equation of motion for 
the thin layer is the same as if the whole foil were accel- 
erated. Finally, we showed that the energy stored in the 
electrostatic field is comparable to the kinetic energy and 
must be taken into account. For applications, the most 
relevant consequences and differences with respect to the 
simplest LS picture arc that the number of "monoener- 
getic" ions is reduced, so that the actual efficiency may 
be quite lower than given by Eq.([T]), and that also light 
ions in a thin layer at the rear surface (e.g., hydrogen 
impurities) may be accelerated by RPA. 
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